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Abstract. Prom covariant transport theory, a significant ~ 10% light quark 
elliptic flow at RHIC implies an elliptic flow of similar magnitude for charm 
quarks, at moderately large px > 2.5 — 3 GeV. At lower transverse momenta, 
charm quark elliptic flow reduces progressively, reminiscent of the mass ordering 
pattern in ideal hydrodynamics. From the quark flows we predict the elliptic flow 
of D mesons at RHIC via quark coalescence. The large parton opacities needed 
to generate the light quark flow also lead to substantial ~ 40 — 50% secondary 
charm production. 



PACS numbers: 12.38.Mh; 24.85.+p; 25.75.Gz; 25.75.-q 

1. Introduction 

Elliptic flow pQ) v 2 = (cos(2<£)), the second Fourier moment of the azimuthal 
momentum distribution, is an important experimental probe that provides information 
about the opacity [21 El El El of the nuclear medium created in noncentral A + A 
reactions. Heavy flavor elliptic flow is especially interesting because it can tell to what 
degree heavy quarks interact and thermalize. Though heavy quarks experience the 
same partonic environment as light partons, their large mass is expected to hamper 
equilibration, at least at RHIC energies ^/snn ~ 100 — 200 GeV. Studying the 
breakdown of equilibrium for the most abundant charm quarks may shed more light 
on the origin of the remarkable thermalization apparent in the light sector [3 [7| [3 E] . 

In principle the experimentally more accessible heavy flavor spectra can 
also constrain dynamical scenarios. However, surprisingly, current (indirect) 
measurements ^U] from Au + Au at RHIC are compatible with both perturbative 
production without charm rescattering and complete charm equilibration . Charm 
hadron, e.g., D meson, elliptic flow will provide a more decisive experimental test |12| . 

Though several schematic scenarios have been explored [131 114) , up to now there 
is no quantitative calculation of charm elliptic flow in heavy ion collisions at RHIC. 
The aim of this study is to make predictions for Au + Au at -^/sjvw = 200 GeV 
based on covariant parton transport theory [21 El CSl CHI 113 113 IE1 EDI That 
approach successfully describes the magnitude and saturation of V2{pt) in the light 
sector, but requires initial parton densities and/or cross sections that are much larger 
than perturbative QCD estimates [3]. Large opacities are also indicated by pion 
interferometry data |18| . Nevertheless, the origin of such opaque conditions is a puzzle. 

One scenario that alleviates the opacity problem is hadronization via quark 
coalescence [21 1221 QUI 1231 123 13 In the coalescence process mesons form 
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from a constituent quark and antiquark, while (anti)baryons from three (anti)quarks. 
Because comoving constituents are strongly favored, elliptic flow can be amplified 22 
03 ED EH 123 EL reducing the opacities needed |S] to explain the data. Utilizing the 
formulas in ^3, we present predictions for the elliptic flow of D and D s mesons at 
RHIC. 



2. Covariant parton transport theory 

We consider here an inelastic extension of the Lorentz-covariant parton transport 
theory in Refs. |2l 131 IT5 1 1161 118j . in which the on-shell parton phase space densities 
{fi(x,p}} evolve with elastic 2 — > 2 and inelastic 2 — > 2 rates as 

P?0/./i,i = J-a£ [[[(h,kh,i^--h,ih,j) \M\> 2 Z k 3 e i \ 2 6 4 (Pi+P2-p 3 -p±) 

jkl 234 

+ S l (x,p 1 ). (1) 

\A4\ 2 is the polarization averaged scattering matrix element squared, the integrals are 
shorthands for J = J d 3 p a / (2E a ), g% is the number of internal degrees of freedom for 

a 

species i, while / a j = fi(x,p a ). The source functions {Si(x,p}} specify the initial 
conditions. 

Though, in principle, could be generalized for bosons and fermions, or inelastic 
3^2 processes ^||2U], no practical algorithm yet exists (for opacities at RHIC) to 
handle the new nonlinearities such extensions introduce. We therefore limit our study 
to quadratic dependence of the collision integral on /. 

We apply Q to a system of massless gluons (g — 16), massless light (u,d) and 
strange quarks/ antiquarks, and charm quarks/ antiquarks (g = 6) with mass M c = 1.2 
GeV. All elastic and inelastic 2-^2 QCD processes were taken into account: gg — > gg, 

gq -> gq, qq -> qq, qq -> qq, qq -> qq; qq -> gg, gg -> qq, and -> q^. The matrix 

elements for massive quarks were taken from j2H]- As in Refs. [21 E], only the most 
divergent parts of the matrix elements were considered, regulated using a Debye mass 
of n D = 0.7 GeV. 

Thus, for all elastic scatterings, da/dt ~ 1/t 2 — ► l/(t — fin) 2 , and for gg <-> qq with 
massless quarks, da/dt ~ ^/(ut) — > l/(u — ^ 2 D ){t — \j? d ). For charm, qq — > cc to first 
approximation can be considered isotropic, while gg *-* cc can be well approximated 
with the angular dependence da/dt ~ 1/[1 — (1 — 4M 2 /s) cos 2 C m] 2 . 

The total cross section for gg — > gg was taken to be constant, neglecting its weak 
logarithmic energy dependence. Its value fixes the total cross sections for all other 
elastic channels: 

G gq^gq ~ (4/9)c r gg^ gg , a qq^qq = (4/9) o'gg^gg ■ (2) 

It also determines the inelastic total cross sections, which on the other hand, decrease 
at high energy. For processes with only massless partons, the energy dependence is 
through the ratio r = njj/s, while for charm, r and R = M 2 / s both play a role |28j : 

_ 2r 1 + r 1 16r 

G 99—qq — ~27 i + 2r r a sg^gg > a qiqi^qjqj — a gg^gg 

2r 



, , 1 + Vl - 4R ,„ , „ „ N VI - 4i? 



(1 + AR + R z ) In \ - (7 + 3i?) 



'gg^gg 



O-gg^ec = ^9(1 - 4i?) 

o-w-cc = 2436(1 - 4H)(1 + 2i?)VT^4l a„ g (3) 
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In order to generate sufficient Ej elliptic flow at RHIC, we take <J gg ^ gg = 10 
mb, about three times the perturbative QCD estimate. As discussed in (jj, this choice 
approximates contributions from 2 <-> 3 inelastic processes |17| . and also relies on the 
amplification of elliptic flow during hadronization via the coalescence process. 

The parton initial conditions for Au + Au at y/s = 200A GeV at RHIC with 
b = 8 fm (w 30% central) were the same as in jj^], except that initial charm 
production was, of course, included. Leading order pQCD minijet three-momentum 
distributions were used (with a if-factor of 2, GRV98LO PDFs, and Q 2 =p?, while 
Q 2 = s for charm). The low-px divergence in the light-parton jet cross sections was 
regulated via a smooth extrapolation below p± < 2 GeV to yield a total parton 
dN(b=0)/dy = 2000 at midrapidity. This choice is motivated by the observed 
dN c h/dy ~ 600 and the expectation that hadronization is dominated by quark 
coalescence [221 123 [Ml El CHI • The transverse density distribution was proportional 
to the binary collision distribution for two Woods-Saxon distributions, therefore 
dN parton (b=8 im)/dy « 500. Perfect r\ = y correlation was assumed. 

The transport solutions were obtained via Molnar's Parton Cascade algorithm 29 
(MPC), which employs the parton subdivision technique |3()j to maintain Lorentz 
covariance and causality. 



3. Results for charm quarks 

Because of the inelastic channels, (QJ is suitable for studying flavor equilibration at 
RHIC (with the limitation that the total parton number is fixed). Results for the 
initial and final rapidity distributions of each parton species are shown in the left 
and right panels of Fig. ^ respectively. The initial condition is dominated by gluons, 
roughly half of which fuse to q — q pairs. This fills in the dip in the light (u, d) 
quark distributions near midrapidity, doubles the light antiquark distributions (u, 
d), enhances strangeness about five-fold, and increases charm by 40 — 50%. Flavor 
changing processes involving only quarks qiq~i — > qjqj influence the flavor chemistry by 
only 10 - 15%. 
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Figure 1. Initial (left panel) and final (right panel) parton rapidity distributions 
in Au + Au at ^/sjvjv = 200 GeV with 6 = 8 fm at RHIC, computed from the 
transport model MPC VM . 



These results show a high degree of light favor and strangeness equilibration at 
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RHIC. The final light antiquark {u, d) and strange quark distributions (s, s) are 
nearly equal, though are still below the light quark (it, d) yields. The baryon number 
distribution, which is one-third of the u + d — (it + d) difference, does not vanish even 
at midrapidity, showing that a baryon free region is not realized at RHIC energies. 

On the other hand, the elliptic flow evolution is driven mainly by elastic 
scatterings. Scatterings off gluons gX — > gX are the most relevant because gluons are 
most abundant and also have larger cross sections J2J- FigureEJi shows the elliptic flow 
of partons at freezeout as a function of pr- For light partons, V2(pr) reaches « 10%, 
in quantitative agreement with j^j . However, instead of saturation between px ~ 1 — 6 
GeV, this calculation shows a turnover and slow decrease of at high pr. The reason 
for the difference is that Ref. |3] considered thermal initial conditions, while here the 
initial spectra are not thermal. The slight difference in Fig. [2K between strange and 
light quark vi is due to the flavor dependence of the initial minijet distributions. 
Because of the larger gluon cross sections (J2J, the gluon V2 exceeds the quark flows. 




p T [GeV] p T [GeV] 



Figure 2. Final parton elliptic flows (left panel) and D meson elliptic flow 
(right panel) as a function of in Au + Au at ^/sjvjv = 200 GeV with 6 = 8 
fm at RHIC, computed using the transport model MPC 1291 . The D meson V2 
was obtained from the parton flows via the quark coalescence formula in 1131 . 
The dashed lines correspond to the two extreme scenarios in [j~3| (see text). For 
hadronization via independent fragmentation, charm hadron vi is approximately 
the same as the charm quark V2- 

At low pt < 2 GeV charm quark elliptic flow is suppressed relative to the V2 of 
light partons, confirming expectations that heavy quark momenta are more difficult 
to randomize. This behavior is reminiscent of the reduction of elliptic flow for heavy 
particles found in ideal hydrodynamics [71 1811 1%1 In contrast, at high pr > 3 
GeV, charm quark V2 is about the same as the light quark V2 ■ At large momenta, the 
heavy-light difference should disappear because the particle mass cannot play a role 
as E/M — > oo. Remarkably, the transport solutions show that most of the difference 
disappears already by pr ~ 2.5 — 3 GeV, i.e., when px/M c ~ a few. 

4. Charm hadron elliptic flow 

From the elliptic flow of partons, one can predict the hadron elliptic flows using a 
suitable hadronization model. Here we are interested in the elliptic flow of the D and 
D s , which will likely be the first charm hadrons measured with sufficient statistics at 
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high p T at RHIC. 

A simple hadronization model is independent fragmentation, in which each parton 
fragments independently to hadrons. In the collinear approximation commonly 
employed, hadron momenta are related to the parton momentum via pr,h = 
zpr,parton, where z is distributed according to the fragmentation function D(z). 
Because for heavy quarks D(z) strongly peaks around z « 1 |32| . charm hadron 
elliptic flow is essentially identical to the charm quark v 2 shown in Fig. [5^,. 

On the other hand, there are strong indications [23 EH IEE] that at intermediate 
2 ^ Pt < 5 GeV hadronization may dominantly occur via quark coalescence [21 1251 
1221 151 125 ) . In the quark coalescence approach, a constituent quark and an antiquark 
that are close in phasespace can combine to form a meson, a(3 — > M, while three 
constituent (anti)quarks can form an (anti)baryon, af3"f — ► B. 

In the simplest (but most successful) variant of the model [221 HI G3) , hadron 
elliptic flow is approximately the sum of constituent flows 

V2,b{x,pt)~ ^ V2,i(x,p T ,i) , v 2 ,m(x,pt) « ^ v 2 ,i{x,p T ,i) , (4) 




with corrections 0({i>2,i} 3 ) that are small in our case. Here J^PT.i = Pt, and the 
hadron momentum is shared roughly in proportion to constituent mass .13;. For 
hadrons composed of u, d, and s quarks, the sharing is approximately equal, while for 
D mesons or the A c , the heavy quark carries most of the momentum. 

Figure Et> shows D and D s elliptic flow as a function of pr given by QJ, for 
constituent mass ratios m Ui< j : m c = 1:5, and m s : m c = 1:3. The results 
are compared to two extreme scenarios |13j : i) zero charm v 2 , i.e., no rescatterings 
at all (lower dashed lines); and ii) v 2 harm (pT) — v 2 9ht (PT) (upper dashed lines), 
which is equivalent to the assumption that D flows the same way as all light mesons 
v^ipr) ~ v 2 (pt) ~ v^ipr)- Neither of the two extremes applies in general. At high 
Pt > 2.5 — 3 GeV, v 2tC ~ v 2 . q and hence the situation agrees with scenario ii), while 
at low pt it is closer to i) because of the suppression in the charm quark v 2 - Thus, D 
meson elliptic flow is predicted to rise smoothly at low pr, and to saturate at w 50% 
higher pt than the pion and kaon v 2 but at the same magnitude. 

The above results ignore fragmentation contributions and therefore are not valid 
at very high pt- For light quarks, the region of validity was estimated 1201 to be 
Pt < 4 — 5 GeV. For charm, the window is smaller because: i) while for light quarks 
coalescence "amplifies" transverse momenta by factor two to three, for charm the 
hadron pt is only 20 — 30% larger than the charm quark pt', ii) fragmentation is 
harder for charm quarks than for lighter partons, therefore charm hadron spectra 
from fragmentation do not fall as steeply at high pt\ and iii) radiative energy loss (jet 
quenching) is expected to be smaller for charm quarks [371 138j than for light partons, 
therefore the fragmentation contribution is less suppressed at high pt- Thus, it is not 
likely that the coalescence results are valid above pt ~ 3 — 4 GeV. 

Finally we emphasize that, despite the success of the simple quark coalescence 
models, severe problems arise when spacetime inhomogeneities or dynamical aspects 
are considered |191 126) . For example, surface- like emission of high-pr particles (as 
seen from covariant transport theory) results in large, spatially nonuniform local 
momentum anisotropics, for which the scaling requires a high degree of lucky 
cancellations to occur |2fi| . 

Furthermore, in the dynamical coalescence approach of Ref. 19 , elliptic flow 
scaling J3J is violated for light quarks. That model combines covariant parton 





i— a,/3,7 



i—a,j3 
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transport theory with the coalescence formula |39| by Gyulassy, Frankel and Remler 
that is applicable to diffuse 4D freezeout distributions in spacetime. Constituents 
without a coalescence partner are assumed to fragment independently. Because 
a significant fraction of constituents has no partner near enough in phasespace to 
coalesce with, elliptic flow is reduced relative to the flow scaling expectation. 

The same problem would be present for charm quarks in that approach, as 
illustrated in Fig. 01 The fraction of light and charm constituents that hadronize 
via coalescence as opposed to fragmentation are shown as a function of pr- Above 
Pt > 3 GeV, less than 25% of charm quarks comes from coalescence in the dynamical 
approach, resulting in a reduction of D meson elliptic flow by almost a half above 
p T » 3.5 GeV. 




1 2 3 4 5 

Pt [GeV] 

Figure 3. Fraction of charm quarks that hadronize via coalescence as opposed 
to fragmentation, as a function pt in Au + Au at ^/sjvjv = 200 GeV with 6 = 8 fm 
at RHIC. Computed in a dynamical coalescence approach 1191 using the transport 
model MPC | 29 | and the Gyulassy-Frankel-Remler coalescence formula 39 . 



5. Conclusions 

This study is the first calculation of charm elliptic flow in Au + Au at ^snn ~ 200 
GeV at RHIC from covariant parton transport theory. Both elastic 2^2 and 
inelastic 2 — > 2 parton interactions were included using MPC [2HI- We show that 
parton opacities needed to generate a w 10% light quark elliptic flow also imply a 
charm quark elliptic flow of similar magnitude, at moderately large pt > 2.5 — 3 
GeV. At lower transverse momenta charm v-i is suppressed, analogously to the mass 
dependence of elliptic flow found in ideal hydrodynamics. Based on hadronization 
via quark coalescence, we predict that the elliptic flows of D and D s saturate at the 
same magnitude « 20% as pion and kaon V2 but at pt ~ 3 GeV, i.e., 50% higher 
than the "meson value" of 2 GeV. In contrast, from hadronization via independent 
fragmentation, charm hadron elliptic flow is only « 10% (same as charm quark V2)- 

We emphasize that the above quark coalescence prediction relies on the 
approximate additivity of elliptic flow, which agrees well with observations in 
the light sector but is quite problematic to preserve in a dynamical coalescence 
approach [23|2E]. In any case, above pt ~ 3 — 4 GeV, fragmentation contributions 
are expected to dominate and reduce D and D s V2 to s» 10%. 
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Finally, the large charm elliptic flow is accompanied by significant flavor 
equilibration at RHIC, with an about five- fold enhancement of strangeness and a 
40 — 50% secondary production of charm. 
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